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Abstract. We review the evidence for cosmic ray acceleration in the superbubble/hot phase of the interstellar medium, 
and discuss the implications for the composition of cosmic rays and the structure and evolution of the interstellar medium 
(ISM). We show that the bulk of the galactic supernovae, their expanding remnants, together with their metal-rich grain 
and gas ejecta, and their cosmic ray accelerating shocks, are all confined within the interiors of hot, low-density su- 
perbubbles, generated by the multiple supernova explosions of massive stars formed in giant OB associations. This 
superbubble/hot phase of the ISM provides throughout the age of the Galaxy a cosmic ray source of essentially con- 
stant metallicity for acceleration by the shocks of many supernovae over time scales of a few Myr, consistent with both 
Be/Fe evolution and ACE observations of ^^Ni/^^Co. We show that the expected metallicity (> 2 times Solar) and filling 
factor (> 50%) of the superbubble/hot phase is high enough that the composition of cosmic rays accelerated from fast, 
supernova grains in these superbubbles is quite consistent with the both Be/Fe and cosmic ray data, while their acceler- 
ation from grains in the well-mixed cooler phases of the ISM is not consistent with observations. We also show that if 
the refractory cosmic ray metals come from the sputtering of fast refractory grains then the accompanying scattering of 
ambient gas by these fast grains can also account for the relative abundance of cosmic ray volatiles. 



INTRODUCTION 

ACE measurements by Wiedenbeck et al. (53) of the 
cosmic ray ^^Ni/^^Co abundance ratio and optical mea- 
surements of the Be/Fe abundance ratios in old stars, e.g. 
Molaro et al. (34) and Boesgaard et al. (4), have shown us 
that the making of cosmic rays is like the making of fine 
wine. Both have to be aged and blended. Cosmic rays 
can not be too young and fresh, because their low ratio 
of K-capture ^^Ni to its daughter ^^Co requires an age of 
> 10^ years before acceleration, and they can not be too 
old and diluted, because the roughly constant ratio of cos- 
mic ray produced Be to supernova produced Fe requires 
the metallicity in the matter from which the cosmic rays 
are accelerated to be roughly constant within a factor of 
~2. This can be achieved if the metallicity is dominated 
by supernova ejecta. 

The cellar where all of this aging and blending hap- 
pens quite naturally is the vast metal-enriched super- 
bubble/hot phase of the interstellar medium (SB/HISM), 
where most supernovae occur. This can be seen from 
the extensive studies of Galactic and extragalactic super- 
novae, their progenitors and their cosmic ray accelerating 
remnants, and of the chemical enrichment and evolution 
of the interstellar medium, that we will review here with 



special emphasis on the metallicity and filling factor of 
the SB/HISM. We also discuss how ACE and other mea- 
surements of the cosmic ray composition can define these 
two properties in the region of cosmic ray acceleration. 

The source of energy for cosmic ray acceleration is 
thought to be shock waves driven by the expansion energy 
of supernova ejecta, e.g. Blandford & Ostriker (3) and 
Axford (1). The power required to maintain the Galac- 
tic cosmic rays is about 10"*^ ergs^', e.g. Lingenfelter 
(23). The Galactic supernova rate is about 3 supernovae 
per century, e.g. van den Bergh & McClure (49), most 
of which (80 to 90%) are core-collapse (Type II and Ib/c) 
supernovae of relatively young (< few 10^ yrs) massive 
O and B stars; the remainder are Type la thermonuclear 
explosions of much older accreting white dwarfs. Thus, 
the average cosmic ray energy needed per supernova is 
about lO^'* ergs, which requires an acceleration efficiency 
of about 10% from the blast wave shocks of the super- 
novae, since supernovae all seem to have similar ejecta ki- 
netic energies of about 10^' ergs, e.g. Woosley & Weaver 
(56) and Nomoto et al. (36). 

The source of the particles that are accelerated as cos- 
mic rays and the site of their acceleration, however, are 
still debated. But if the energy comes from supernova 
shocks, the site and source of the particles clearly must 



be in the material through which the shocks pass. New 
clues to origin of the particles come from recent measure- 
ments, e.g. Molaro et al. (34) and Boesgaard et al. (4), of 
Be/Fe abundances in old halo stars that show that the ra- 
tio of cosmic ray spallation produced Be relative to core- 
collapse supernova-produced Fe has remained roughly 
constant throughout the evolution of the Galaxy. This 
constancy requires (Ramaty, Kozlovsky & Lingenfelter 
(38); Ramaty, Lingenfelter & Kozlovsky (40)) that the 
cosmic rays be accelerated out of matter that is only 
partially diluted by mixing with the interstellar medium 
(ISM) so that it is still sufficiently enriched in supernova- 
synthesized metals that its metallicity did not changed by 
more than ~2 over galactic evolution. Detailed calcu- 
lations by Ramaty et al. (42) and Ramaty, Lingenfelter 
& Kozlovsky (41) of the production and evolution of the 
Galactic Be/Fe ratio clearly show that the bulk of the cos- 
mic rays can not be accelerated from the well-mixed ISM, 
as has been recently assumed, e.g. Meyer et al. (33) and 
Ellison etal. (11). 

Just such a metal-enriched environment is where most 
supernovae do in fact occur and through which most of 
the cosmic ray accelerating, supernova shocks propagate. 
As we have discussed in Higdon, Lingenfelter & Ramaty 
(17,18), extensive observations show that the bulk of the 
core-collapse supernova progenitors are formed in OB as- 
sociations in giant molecular clouds and that the com- 
bined winds and supernova ejecta of these stars form hot, 
low density superbubbles, that reach dimensions of sev- 
eral hundred pc and last for tens of Myr During this time 
the bulk of the supernova ejecta, the supernova shocks, 
and their cosmic ray acceleration are all confined within 
the superbubble/hot phase of the ISM. 



SUPERNOVAE & THE 
SUPERBUBBLE/HOT ISM 

Core-collapse (Type II and Ib/c) supernovae are highly 
correlated in space and time, e.g. McCray & Snow (28). 
Such supernovae thus create giant cavities, or superbub- 
bles, in the interstellar medium rather than many smaller, 
isolated bubbles, e.g. Mac Low & McCray (26) and 
Tomisaka (48). This is expected because 1) the massive O 
and B star supernova progenitors (>8 Mq; e.g. Woosley 
& Weaver (56); Nomoto et al. (36)) are not distributed 
uniformly in interstellar space, but tend to form clusters, 
since the majority of these massive stars are born in the 
most massive (> 10^ M,:,) molecular clouds in gravi- 
tationally unbound OB associations, while less massive 
clouds are destroyed by the intense UV irradiation with 
the birth of their first O star (McKee & Williams (30)); 
and 2) these stars are short-lived and slow moving; the 



progenitors of core-collapse supernovae have main se- 
quence lifetimes of ~ 3 to 35 Myr and OB stars in as- 
sociations have dispersion velocities of only ~4 km s^' 
(Blaauw (2)), so they do not travel too far (^^120 pc in 30 
Myr) from their birthplaces before they die in supernova 
explosions. Consequently, the combined effect of these 
clustered supernova explosions is to create superbubbles, 
which expand and merge to form the hot (>10^ K), ten- 
uous (<10^^ cm^^) phase of the ISM with an average 
filling factor of ~ 50%, or more (e.g. Yorke (57); Spitzer 
(45); McKee (29); Rosen & Bregman (43); Korpi et al. 
(20)) of the Galactic disk (with a scale height ^^100 pc) 
and essentially all of the corona/halo (with scale height 
of ^3 kpc). We discuss the spatial distribution of the 
SB/HISM filling factor in more detail below. 

An analysis of the surface brightness distribution of 
the remnants of historical supernovae in our Galaxy by 
(Higdon & Lingenfelter (16) has shown that 85 ±10%, 
of the observed Galactic supernovae occured in the su- 
perbubble hot phase of the ISM. This is quite consistent 
with more extensive observations of supernovae in other 
late type galaxies. As discussed in detail in Higdon et al. 
(17), the combined observations of van Dyk et al. (51) 
and Kennicutt, Edgar, & Hodge (19) show that the great 
majority, ~ 90±10%, of the core-collapse supernovae in 
late type galaxies also occur within superbubbles, and be- 
cause of the large filling factor of the superbubble, hot- 
phase of the ISM, half, or more, of the Type la should 
also occur within the superbubbles just by chance. Thus, 
since core-collapse supernovae account for 80 to 90% of 
all supernovae in our galaxy and Type la make up the 
remainder, roughly 80% of all supernovae occur in the 
superbubble hot phase of the ISM, and the bulk of the 
cosmic rays accelerated by their shocks are also produced 
there. 

The observed concentration of supernovae in the 
SB/HISM and the subsequent cosmic ray acceleration in 
this hot (> 10^ K) phase also argues strongly against a 
first-ionization-potential (FIP) injection bias for cosmic 
ray enrichment which requires warm partially ionized 
gas, and does not work in the nearly fully ionized gas 
of the SB/HISM. However, the SB/HISM environment is 
quite consistent with cosmic ray source particle injection 
from the sputtering of the high velocity refractory grains 
formed in supernovae, see Lingenfelter, Ramaty & Ko- 
zlovsky (25) and Lingenfelter & Ramaty (24). Such a 
volatility bias for cosmic ray refractory metal injection 
from sputtering of supernova grains was first suggested 
by Cesarsky & Bibring (7) and it was also recently pro- 
posed by Meyer et al. (33) for shock-accelerated, ice- 
stripped, refractory cores of grains in the warm ISM. 

The occurrence of most supernovae in the SB/HISM 
and cosmic ray acceleration in that hot phase also argues 
strongly against the mass/charge (A/Q)-dependent accel- 



eration model of Ellison et al. (11) for the volatile ele- 
ments, because this too requires warm partially ionized 
gas, and Ellison & Meyer (12) argue it does not work 
in the highly ionized gas of the SB/HISM. As we have 
shown in Lingenfelter & Ramaty (24) and discuss below, 
however, a mass dependent injection of volatiles appears 
to result quite naturally from the scattering of the ambi- 
ent gas atoms in direct collisions with fast grain atoms 
that must accompany the sputtering of the grains. 



METALLICITY OF THE 
SUPERBUBBLE/HOT ISM 

The bulk of the metals (elements with Z>5) in the 
Galaxy have been produced by supernovae and ejected 
into the ISM. The relative abundances of most elements 
have remained relatively constant (e.g. Timmes, Woosley 
& Weaver (47)), because they simply reflect the IMF- 
averaged supernova yields which do not depend strongly 
on the interstellar metallicity, e.g. Woosley & Weaver 
(56). The averaged relative abundances of the present 
ISM, e.g. Savage & Sembach (44), the older (4.5 Gyr) 
Solar system material, e.g. Grevesse, Noels & Sauval 
(15), and the IMF-averaged fresh supernova ejecta, e.g. 
Lingenfelter et al. (25), are all within about 10% of 
one another But their overall abundance in the ISM (i.e. 
the interstellar metallicity) has grown steadily over time 
with the accumulation of fresh supernova ejecta continu- 
ously injected and mixed into the ISM. The time scale for 
thorough mixing is generally thought (e.g. McWilliam 
(31); Thomas, Greggio & Bender (46)) to be on the or- 
der of 30 to 100 Myr This is comparable to the typical 
mean life of the SB/HISM reservoir into which the bulk 
of the supernova ejecta with a metallicity 10 times Solar 
(e.g. Woosley & Weaver (56)) is injected and in which the 
bulk of the mixing is expected to occur. Thus, we would 
expect significant variations in the metallicity, but not in 
the abundances of most elements relative to one another, 
within the SB/HISM as a function of the age of individual 
superbubbles and their generating OB associations. 

The average, or equillibrium, metallicity of the 
SB/HISM is not known, but the supernova ejecta appear 
to be able to provide sufficient metals to produce a metal- 
licity >2 times Solar, as is required (Ramaty et al. 2000b) 
for the cosmic ray source from the constancy of the Be/Fe 
abindances in old stars. This can be seen from a sim- 
ple comparison of the total mass of the SB/HISM and 
the mass of supernova ejecta produced during the mean 
life of the SB/HISM. The total mass of the SB/HISM 
is ^ 10^ M(V), assuming an average, e.g. Yorke (57) 
and Spitzer (45), SB/HISM density of - 10"^ H/cm-\ a 
SB/HISM scale height of ~ 3 kpc and an effective Galac- 



tic radius of 15 kpc. Taking a nominal SB/HISM 
mean life, or mixing time, f 100 Myr, the required 
SB/HISM input is - 1 M0/yr(f/lOOMyr). The present 
Galactic SNII/Ibc rate of about 1 SN every 40 yr, pro- 
ducing an IMF-averaged ejecta mass of 18 M©, gives a 
Galactic SNII/lbc ejecta input of ^ 0.45 M0/yr with a 
metallicity, zsn of 10 times Solar. If all of the remain- 
ing SB/HISM mass comes from evaporated clouds and 
swept up gas in the well-mixed ISM with Solar metal- 
licity zq of 1, then the averaged SB/HISM metallicity, 
zmsM ~ ['^OMsn + '^{Mhism-Msn)]/Mhism, or zhism ~ 
1 +9Msn/Mhism 1 H-4(f/100Myr). Thus the SB/HISM 
metallicity zhism > 2 times Solar for any SB/HISM mean 
life, or mixing time, f > 25 Myr, consistent with the esti- 
mated values, e.g. McWilliam (31) and Thomas, Greggio 
& Bender (46). 

Such a mixing time, or mean age, of metals from su- 
pernova ejecta in these superbubbles is more than a cou- 
ple orders of magnitude longer than the minimum age 
(<100 kyr) of cosmic ray source metals required by the 
ACE observations of Wiedenbeck et al. (53), showing 
that the bulk of the ^^Ni had decayed (with a 110 kyr 
mean life) in the cosmic ray source material prior to ac- 
celeration. 

Observational evidence of such supernova ejecta en- 
riched superbubble metallicity may be found in the x-ray 
emission from the interiors of giant HII regions in the 
Large Magellanic Cloud (LMC), thought to be superbub- 
bles powered by supernovae. The observed x-ray lumi- 
nosities of these bubbles, which should scale directly with 
metallicity, are an order of magnitude higher than would 
be expected (Chu & Mac Low (8)) if they had a typical 
LMC metallicity of only 1/3 Solar. Thus, we suggest that 
the x-ray observations do in fact imply an average metal- 
licity of roughly 3 times Solar in these superbubbles. 

Such metallicities are larger than that calculated from 
the simple analytic superbubble model of Mac Low & 
McCray (26), which assumes conductive heating and 
evaporation of swept-up ISM as the primary source of 
superbubble gas and predicts an averaged metallicity of 
only 1 . 1 for a ~ 50 Myr old superbubble of ~ 700 pc 
radius. But this model was based on several assumptions 
that greatly reduce the superbubble metallicity. First, 
the model neglected the interstellar magnetic fields that 
would greatly supress the conductive heating normal to 
the field lines, overestimating the ISM input, and also 
provide additional external confining pressure, underesti- 
mating the supernova input required to generate the bub- 
ble, e.g. Tomisaka (48). In addition, the model assumed 
the unit density ISM extended to heights much larger than 
the 700 pc bubble radius, instead of the measured scale 
height of < 200 pc, which would greatly reduce the as- 
sumed ISM input, and moreover would allow the blow- 
out of the superbubble into the Galactic halo, also greatly 



underestimating the supernova input required to generate 
the bubble. As a result of these underestimates of the 
required supernova power, the model was able to gen- 
erate a ~ 50 Myr old superbubble of ~ 700 pc radius, 
with a very low effective rate of SNII/Ibc supernovae of 
only ^ 2 SN/Myr kpc^ in the Galactic plane. This as- 
sumed rate is only 3% of the estimated local SNII/Ibc 
rate ~ 70 SN/Myr kpc^ (as we show below), and thus 
the model underestimates the supernova ejecta input by a 
factor of 35 ! When an appropriate supernova rate is used 
and even minimal effects of the magnetic fields and gas 
scale height are considered, superbubble metallicities of 
more that 2 times Solar would be expected. 

A large fraction of the C, O and refractory metals in 
this ejecta may be in graphite and oxide grains, since 
in the core-collapse supernova 1987A roughly 0.2 Mq 
of this material condensed out of the cooling, expand- 
ing ejecta as high velocity 2500 km/s) grains within 
2 years after the explosion, see Kozasa, Hasegawa & 
Nomoto (21), and as much as 1 Mq could be expected, 
see Dwek (9), to condense before the ejecta is reheated 
and slowed by the reverse shock and the grains with a 
much smaller charge to mass ratio begin to move sepa- 
rately from the ejecta plasma. In fact, Dwek (9,10) sug- 
gests that supernova ejecta are the major source of refrac- 
tory grains in the Galaxy and interactions with supernova 
shocks are the major cause of their destruction. 

Thus, supernova ejecta and winds can be expected to 
dominate the metallicity and grains within the SB/HISM, 
where the bulk of supernova shock waves are dissi- 
pated and the bulk of cosmic rays should be acceler- 
ated. These supernova grains should therefore be the ma- 
jor injection source required for the cosmic ray metals, 
because of their high initial velocity (Lingenfelter et al. 
(25)) and possible subsequent acceleration (Elhson et al. 
(11)). Moreover, because the metallicity of the supernova 
ejecta is essentially independent of progenitor metallic- 
ity (Woosley & Weaver (56)), the SB/HISM can provide 
the essentially constant source of cosmic ray metals re- 
quired by Be/Fe observations. Therefore, we would ex- 
pect that throughout the age of the Galaxy, the bulk of 
the core-collapse supernovae occur in the metal enriched 
SB/HISM, and the blast wave shocks of their remnants 
accelerate the bulk of the Galactic cosmic rays out of the 
enriched gas and dust in the SB/HISM. 



FILLING FACTOR OF THE 
SUPERBUBBLE/HOT ISM 

The hot (- 10*^ K), tenuous (~ 10"^ H/cm^) phase of 
the interstellar medium is powered primarily by Galac- 
tic supernovae and formed through the merger of super- 



bubbles, generated by the clustered supernovae in OB 
associations. This can be seen energeticly from a com- 
parison of the power required to maintain the pressure 
in the SB/HISM and that provided by Galactic super- 
novae, which suggests that the filling factor, i.e. the frac- 
tional volume, of the SB/HISM should be large. The to- 
tal energy in SB/HISM is ~ 3 x IQ^^ fnisM ergs, assum- 
ing a SB/HISM filling factor, fyisM, a SB/HISM pres- 
sure of 3 X 10^'^ erg/cm^, a Galactic radius of 15 
kpc and a scale height of 3 kpc. For a SB/HISM mean 
life t of 100 Myr, the power required to maintain the 
SB/HISM is - 3 X 10'^^ fmsM ergs/yr(f/100Myr). The 
Galactic SNII/Ibc rate of ^ 1 SN/40 yr with an average 
ejecta energy ^10^' ergs/SN, gives a Galactic SNII/Ibc 
power of 2.5 X 10"*^ ergs/yr Thus even with signifi- 
cant (>50%) energy losses SNII/Ibc could completely fill 
(hiSM = 1) the Galaxy with the SB/HISM in r > 25 Myr. 

The overall Galactic average value of the filling factor 
of SB/HISM is, in fact, generally taken to be ~ 50%, or 
more, depending on the assumed Galactic scale height, 
e.g. Yorke (57), Spitzer (45), McKee (29). Because of 
the strong dependence of the SB/HISM filling factor on 
local supernova rates, it is thought to be high ~ 90% in 
the inner Galaxy (i.e. within the Solar radius of 8.5 kpc) 
where most Galactic supernovae occur, as well as in the 
Galactic halo where the superbubbles blow-out, and low 
< 50% in the outer Galaxy beyond the Solar radius where 
few supernovae occur 

A more quatitative estimate of the dependence of the 
SB/HISM fining factor on Galactic radius (see Table 1) 
can be made from the radial dependence of the Galactic 
supernova rate and the calculated filling factor versus the 
supernova rate. 

The dependence of the SB/HISM filling factor on local 
supernova rates has been quantified by recent calculations 
of 2D hydrodynamics by Rosen & Bregman (43), and 3D 
magnetohydrodynamics by Korpi et al. (20). These cal- 
culations determined the filling factors of all phases of the 
ISM as a function of height z above the Galactic plane for 
a range of supernova rates. Generally these calculations 
suggest that the SB/HISM filling factor is lowest at the 
Galactic plane where the superbubble expansion is most 
constrained by the warm and cold phases of the ISM, and 
increases to ~ 100% in the halo at large distances above 
the plane. For the purposes of cosmic ray acceleration 
and composition, what is important is the height averaged 
fining factor of the SB/HISM for | z |< 300 pc, which 
is the range of heights where most of the supernovae 
occur. For assumed supernova rates (adjusted to 10^' 
ergs/SN) of 5, 20, 40 and 80 SN/Myr kpc^ in the Galactic 
plane these calculations give height averaged (| z |<300 
pc) SB/HISM filling factors of ~ 0.1, 0.4, 0.6 and 0.9 
respectively. This suggests that at low supernova rates 
(< 40 SN/Myr kpc^) the SB/HISM filling factors within 



Table 1. RADIAL DEPENDENCE OF GALACTIC SUPERNOVA 
RATE & EXPECTED SB/HISM FILLING FACTOR 



Galactic 


MoleCloud 


OBAssoc 


Supernova 


SB/HISM 


Radius 


Density 


Density 


Rate * 


Filling Factor' 


kpc 


MO/pc^ 


N/kpc- 


SN/kpc-Myr 


z|<300pc 


1 


1 


0.3 


35 


~0.4 


4 


7 


2.3 


250 


~0.9 


6 


6 


1.6 


210 


~0.9 


8 


2.5 


0.6 


90 


~0.9 


10 


1.5 


0.4 


50 


~0.5 


15 


0.4 




12 


~0.1 



* Galactic SN rate of 3 SN/lOOyr normalized to surface density distribution of 
molecular clouds from Williams & McKee (54) and OB associations from 
McKee & Williams (30) 

" Expected filling factor based on the SN rate from the hydrodynamic 
calculations by Rosen & Bregman (43) and Korpi et al. (20). 



z |<300 pc scale roughly linearly with the supernova 
power, while at higher supernova rates (> 80 SN/Myr 
kpc2) the SB/HISM filling factors within | z |<300 pc 
reach a maximum value of ^ 90%. 

The Galactic radial dependence of the supernova rate 
can be estimated by normalizing the Galactic SN rate of 
3 SN/lOOyr to the radial dependence of the surface den- 
sity of either molecular clouds from Williams & McKee 
(54) or OB associations from McKee & Williams (30), 
which are proportional to one another, as we see in Table 
1 . Such a normaUzation gives a local supernova rate at the 
Solar distance (8.5 kpc) of about 80 SN/Myr kpc^ and a 
peak rate at about 4 kpc of 250 SN/Myr kpc'. From the 
calculated dependence of the filling factor on supernova 
rate, we thus estimate the Galactic radial dependence of 
the SB/HISM filHng factor within | z |<300pc, as shown 
in Table 1. We see that the SB/HISM is expected to fill 
most (-- 90%) of the ISM within | z | <300 pc from some- 
where inside of 4 kpc out to roughly the Solar distance 
of 8.5 kpc, decreasing thereafter with Galactic radius to 
~ 50% at 10 kpc and ^ 10% at 12 kpc where the super- 
nova rate is very low. As we show below, a SB/HISM 
filling factor of > 50% can provide a cosmic ray injec- 
tion composition in the SB/HISM that is consistent with 
current estimates of the required cosmic ray source com- 
position. We note that one recent estimate by Ferriere 
(14) of the radial dependence of the SB/HISM filHng fac- 
tor gives only 20% locally, but this is for a very low local 
supernova rate from a very steep assumed radial depen- 
dence that is not consistent with the molecular cloud and 
OB association observations. 

A local SB/HISM filHng factor of - 90% would ap- 
pear to be quite consistent with observations within the 
local kpc, see Blaauw (2) Fig. 8, which show that the 



Sun presently lies inside the ~ 500 pc radius superbub- 
ble produced by the ~ 30 Myr Cas-Tau OB association, 
e.g. Olano (37). This local superbubble is defined in the 
Galactic plane by a ring of young OB associations know 
as Gould's Belt which have formed from the ring of cool- 
ing gas swept up by the superbubble. The Cas-Tau as- 
sociation inturn is part of a larger 1 kpc radius) ring 
of OB associations, including Cam-1, Aur-1, Gem-1 and 
Mon-2, formed by an older, now vanished OB associa- 
tion. 



COSMIC RAY ACCELERATION IN 
SUPERBUBBLE/HOT ISM 

These hot, low density superbubbles are the hot phase 
of the ISM, where shock acceleration of cosmic rays is 
expected, e.g. Axford (1), to be "most effective", be- 
cause the energy losses of the accelerated particles are 
greatly reduced and the supernova shocks do not suffer 
major radiative losses, as they would in a denser medium. 
The rapid radiative loss of supernova remnant energy in 
the average ISM sets in at a radius of ^^20 pc, while the 
undiminished shock energy of nonradiative remnants in 
the superbubble hot phase expand out to radii of ~200 
pc. At full shock energy, supernovae in the low density 
SB/HISM expand to ^10^ times the volume of those in 
the average ISM. Thus, the supernova shocks in low den- 
sity, but metal enriched SB/HISM process a comparable 
masses of gas and for z > 2 at least twice the metals as 
those in the average ISM, contrary to the estimate of EI- 
Hson & Meyer (12). 

Also, since the energy of supernova shocks in the 
SB/HISM, unlike that of shocks in the denser ISM, is 



not dissipated by radiation losses before the shocks slow 
to sound speed, cosmic rays are accelerated in SB/HISM 
primarily by low Mach number shocks. Such low Mach 
number (e.g. <4) shocks can produce, e.g. Axford (1), 
the power-law index of ~2.3 required for the cosmic ray 
source spectrum, while the lower spectral indices (^2) 
produced by high Mach number shocks in the denser ISM 
are not consistent with the required source value. 

The observed concentration of supernovae in the su- 
perbubble hot phase and the much higher acceleration 
efficiency expected there clearly show that the bulk of 
the cosmic rays must be accelerated in the SB/HISM. 
Such an acceleration site also argues strongly against a 
first-ionization-potential (FIP) injection bias, e.g. Meyer 
(32), which requires warm partially ionized gas, not the 
highly ionized gas of the hot phase. Acceleration in 
the SB/HISM further argues against a mass/charge (A/Q) 
dependent acceleration model for the volatile elements, 
which Ellison & Meyer (12) argue does not work in 
highly ionized hot gas. As we have shown in Lingenfelter 
et al. (25) and Lingenfelter & Ramaty (24) and discuss 
further below, however, sputtering and scattering of hot 
gas by high velocity refractory grains from supernovae 
in the SB/HISM can provide a self-consistent cosmic ray 
injection source for both refractory and volatile elements. 

The transient acceleration of low energy (<100 
MeV/nucleon) cosmic rays (LECRs) in superbubbles has 
also been suggested, e.g. Bykov (5), as an alternative 
source of Be production in the Galaxy. To account for 
the measured Be/Fe evolution solely by LECRs, how- 
ever, would require (Ramaty et al. (41)) that there be as 
much or more energy in the LECRs as there is the rela- 
tivistic cosmic rays. Bykov (5) suggests that such LECRs 
might be accelerated in supernova shocks during the early 
(<3 Myr) stages of superbubble formation and that these 
LECRs are later further accelerated to relativistic cosmic 
ray energies by the ensemble of supernova shocks as the 
superbubble fully develops, e.g. Bykov & Fleishman (6). 
But since the energy in such LECRs persists for only a 
small fraction (<10%) of the age (~50 Myr) of the su- 
perbubble and then more energy is added as the LECRs 
become relativistic cosmic rays which persist for most of 
the age of the superbubble, such a model can not produce 
a time averaged LECR energy comparable to that of the 
relativistic cosmic rays. Even if the LECRs were not fur- 
ther accelerated to relativistic energies, comparable total 
energy densities in LECRs and relativistic cosmic rays 
would require that roughly half of the supernovae in su- 
perbubbles accelerate LECRS, but <5% of the superbub- 
ble supernovae occur during the first few Myr of super- 
bubble growth when condition favorable to LECR accel- 
eration might be expected (Bykov (5)). 



EXPECTED ABUNDANCES OF 
REFRACTORY COSMIC RAYS 

We have shown in Lingenfelter et al. (25), Higdon 
et al. (17) and Lingenfelter & Ramaty (24) that the ob- 
served enrichment of the cosmic ray refractory elements 
can be produced by the preferential acceleration in the 
SB/HISM of suprathermal ions sputtered off high veloc- 
ity (few 1000 km s^') refractory grains, which formed as 
condensates in the expanding ejecta of supernovae, e.g. 
Kozasa et al. (21) and Dwek (9). The measured (Naya et 
al. (35)) broad width (5.4±1.4 keV) of the Galactic 1.809 
MeV line from the decay of long-lived (1.0x10^ yr mean 
life) ^^Al, most likely produced in Type II supernovae, 
e.g. Woosley & Weaver (56), clearly suggests that re- 
fractory grains, containing most of the live Galactic ^^Al, 
are still moving at velocities of ^^450 km s^ some 10^ 
yrs after their formation, and that the bulk of the grains 
are in low density superbubbles because the grains would 
have been stopped much earlier in the much denser aver- 
age ISM. We also showed that only a very small fraction 
(-10-4) of the grains formed in a typical supernova need 
be accelerated to account for the average injection of cos- 
mic ray metals. 

The similarity of the cosmic ray source and solar abun- 
dance ratios of refractory elements, mainly Mg, Al, Si, 
Ca relative to Fe, simply reflects the fact that supernovae 
are the primary source of these elements, e.g. Timmes 
et al. (47), and that the SB/HISM filUng factor is large 
where cosmic rays are accelerated, so that the bulk of the 
Fe grains from the SNIa also contribute to the high veloc- 
ity grain population in the SB/HISM. In particular, since 
the Si, Mg, Al, and other refractory elements are primar- 
ily produced in core-collapse SNII/Ibc, while only about 
half of the Fe is made in them and the other half is made 
in thermonuclear SNIa, a SB/HISM filling factor of — 
90% leads to differences of only ~ 5% between the Si/Fe 
ratio in SB/HISM and the average Galactic production ra- 
tio, which determines that in the well-mixed ISM. This is 
well within the present uncertainties in the inferred cos- 
mic ray source ratios shown in Table 2, where we see 
from a much more detailed estimate in Lingenfelter & 
Ramaty (24) that the injection abundances expected for 
cosmic ray acceleration predominantly in the SB/HISM 
is consistent with the present cosmic ray source ratios of 
Engelmann et al. (13) even for an assumed SB/HISM fill- 
ing factor of only 50%. Similar small differences < 10% 
in relative abundances from a SB/HISM filling factor of 
^ 90% would be expected for those s-process elements 
which appear to come primarily from the winds of less 
massive stars. 

The estimated mean refractory abundances in super- 
nova grains (Table 2) are based on the calculations by 



Table 2. COSMIC RAY INJECTION ABUNDANCE RATIOS IN % 



ISMGrains ISMCores * SNGrains"'" SBGrains** CRInject* CRSource^' Solar' 



C/Fe 


690 


-?- 


210-510 






422±14 


1122±139 


0/Fe 


1400 


400 


320-520 


460-690 


455-665 


522±11 


2344±414 


Mg/Fe 


115 


110 


50-150 


90-190 


90-185 


103±3 


120±4 


Al/Fe 


10 


10 


5-16 


8-20 


8-20 


7.7±1.5 


9.8±0.3 


Si/Fe 


105 


65 


110-170 


105-185 


100-175 


99±2 


115±4 


Ca/Fe 


6 


6 


4-8 


5-9 


5-9 


6.0±0.9 


7.1 ±0.2 


Ni/Fe 


6 


6 


6-14 


6-9 


6-9 


5.6±0.2 


5.6±0.2 



* ISMGrains and ISMCores - HST interstellar depletion determined abundance from Savage & Sembach (44). 
^ SNGrains - Range of IMF averaged supernova ejecta mixes weighted with relative SNII:SNIb:SNIa rates of 
67-75%:13-15%:20-10% from van den Berg & Tammann (50) and van den Berg & McClure (49), except for O; for 
the SNII and SNIb contributions, refractory O is assumed to be bound in MgSiOs, Fe304, AI2O3, CaO and NiO, and 
for (the very small) SNIa contribution, all the produced O is assumed bound to Fe. 

** SBGrains - Modified SNGrains for 85% of SNII and SNIb and 50% of SNIa in superbubbles plus ISM refractory 
grain ISMCores for a mean superbubble metallicity range of 2-5 times that of ISM, as discussed in the text. 

CRInject - Galactic supernova averaged grain abundances for cosmic ray injection, taking a mix of SBGrain 
abundances for supernova acceleration in superbubbles and ISMCore grain abundances (without any supernova 
emichment) for supernova acceleration outside the superbubbles, weighted by the relative swept-up metal masses and 
supernova rates, as discussed in the text. 

^ CRSource - elemental abundances from Engelmann et al. (13). 

" Solar system - elemental abundances from Grevesse, Noels & Sauval (15). 



Woosley & Weaver (56), Woosley, Langer & Weaver (55) 
and Nomoto et al. (36) of supernova ejecta abundances 
for Types 11, lb and la, averaged over the initial mass 
function and supernova rates of van den Berg & Tam- 
mann (50) and van den Bergh & McClure (49), except 
that we assume the grain O abundance is limited to that 
bound in AI2O3, MgSiOa, Fe304, CaO and NiO. We 
also show for comparison, the refractory abundances in 
the typical, older icy interstellar grains (ISMGrains) and 
their refractory cores (ISMCores) recently determined by 
HST observations, see Savage & Sembach (44). Here 
we see that the Si/Fe of 65% in refractory cores of ISM 
grains, which Meyer et al. (33) proposed as the cosmic 
ray source, is not consistent with the required cosmic ray 
source value of 99±2%. 



EXPECTED ABUNDANCES OF 
VOLATILE COSMIC RAYS 

In addition to the sputtering of refractory ions, the in- 
teractions of the high velocity, supernova grains can also 
provide a simultaneous, self consistent cosmic ray injec- 
tion source of H, He and other volatiles. Cesarsky & 
Bibring (7) suggested that high velocity grains may tem- 
porarily pick up by implantation volatile atoms from the 
gas through which they pass, and their subsequent sput- 
tering could provide a source of less enriched suprather- 
mal volatiles. We suggest a much more direct injection 



process for the volatiles. Since direct collisions of fast 
grains with ambient gas atoms and ions are thought to be 
the primary means of grain momentum loss, e.g. Elli- 
son et al. (11) §2.3, we would expect that the supernova 
grains should simply scatter ambient H, He and other 
volatile atoms to the same suprathermal injection veloci- 
ties as the grains and their sputtered refractory products. 
Such a process would, in fact, directly account for the 
measured cosmic ray abundance ratio by number of the 
refractory (including C and "bound" O) to volatile ele- 
ments, i.e. (C,0,Mg,Al,Si,Fe,etc)/(H,He,etc) = 0.010 of 
Engelmann et al. (13), since Ellison et al. (11 §2.4) as- 
sume that roughly 0.5%- 1% of grain collisions with am- 
bient gas atoms, predominantly scattering volatile atoms, 
result in the sputtering of a refractory atom from the grain 
surface, all of which come off with essentially the same 
injection velocity. Moreover, because the geometric scat- 
tering cross section increases with mass to the 2/3 power, 
such scattering should also lead to a mass-dependent en- 
richment of heavier volatiles with respect to H, as is ob- 
served in the cosmic rays, e.g. Meyer et al. (33), and 
which Ellison & Meyer (12) argue can not be accounted 
for by an A/Z dependent acceleration bias in the hot ISM. 

The composition of the grain-scattered suprathermal 
volatiles can be further enriched by the fact that most 
of the supernova shocks will be interacting with grains 
and gas in the supernova-ejecta and progenitor-wind en- 
riched superbubbles. Since the ^^Ne/^'^Ne ratio in the 
Wolf Rayet winds of massive, supernova progenitors may 
exceed the solar system value by more than two or- 



ders of magnitude, e.g. Maeder & Meynet (27), grain- 
scattering of such wind enriched could account for the 
high ^^Ne/^'^Ne observed in the cosmic rays, e.g. Leske 
et al. (22). The existence of such a Wolf Rayet signature 
in the cosmic rays also provides further evidence for the 
acceleration of cosmic rays in the superbubble hot phase 
where the bulk of the massive Wolf Rayet, supernova pro- 
genitors are also confined. 

This work was supported by NASA ATP and ACE/GI 
Programs. 
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